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In a recent paper we used chiral perturbation theory to determine the phase diagram and pion 
spectrum for Wilson and twisted-mass fermions at nonzero lattice spacing with nondegenerate up 
and down quarks. Here we extend this work to include the effects of electromagnetism, so that it is 
applicable to recent simulations incorporating all sources of isospin breaking. For Wilson fermions, 
we find that the phase diagram is unaffected by the inclusion of electromagnetism—the only effect is 
to raise the charged pion masses. For maximally twisted fermions, we previously took the twist and 
isospin-breaking directions to be different, in order that the fermion determinant is real and positive. 

However, this is incompatible with electromagnetic gauge invariance, and so here we take the twist 
to be in the isospin-breaking direction, following the RM123 collaboration. We map out the phase 
diagram in this case, which has not previously been studied. The results differ from those obtained 
with different twist and isospin directions. One practical issue when including electromagnetism is 
that the critical masses for up and down quarks differ. We show that one of the criteria suggested 
to determine these critical masses does not work, and propose an alternative. 


I. INTRODUCTION 

The phase diagram of lattice QCD (LQCD) can con¬ 
tain unphysical transitions and unwanted phases due to 
discretization effects. A well known example is the Aoki 
phase that can be present with Wilson-like fermions mQ 
Unphysical phases occur when the effects of physical light 
quark masses are comparable to those induced by dis¬ 
cretization, specifically m ^ a^Ag^j^, with a the lattice 
spacing. This can be shown by extending chiral pertur¬ 
bation theory (XPT) to include the effects of discretiza¬ 
tion [2]. Understanding the phase structure is necessary 
so that LQCD simulations can avoid working close to un¬ 
physical phases, so as to avoid distortion of results and 
critical slowing down. 

Recently, we extended the analysis of the phase dia¬ 
gram to the case of nondegenerate up and down quarks 
for Wilson-like and twisted-mass fermions [3]. This was 
prompted by the recent incc^oration of mass splittings 
into simulations of LQCDj^ We found a fairly com¬ 
plicated phase structure, in which, for example, the 
Aoki phase was continuously connected to Dashen’s CP- 
violating phase [a [7]. 

A drawback of our analysis was that it did not in¬ 
clude the other major source of isospin breaking in QCD, 
namely electromagnetism. For most hadron properties, 
electromagnetic effects are comparable to those of the 
mass nondegeneracy Cq = {mu — vn^ j2. For example, in 
the neutron-proton mass difference these two effects lead 
to contributions of approximately —1 MeV and 2.5 MeV, 
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^ “Wilson-like” refers to both unimproved and improved versions 
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respectively]^ Furthermore, the recent LQCD simula¬ 
tions alluded to above have included both mass nonde¬ 
generacy and electromagnetism. Thus, to be directly ap¬ 
plicable to such simulations, we must extend our analysis 
to include electromagnetism. This is the purpose of the 
present note. 

We work in Wilson or twisted-mass XPT (both of 
which we refer to as WXPT for the sake of brevity) us¬ 
ing a power-counting to be explained in Sec. [Ilj At the 
order we work, it turns out that the inclusion of elec¬ 
tromagnetism can be accomplished in most cases simply 
by shifting low-energy coefficients (LECs) in the results 
without electromagnetism. Thus we can take over many 
results from Ref. [3] without further work. 

One new issue concerns the simultaneous inclusion 
of electromagnetism and quark nondegeneracy with 
twisted-mass fermions. The approach we used in the ab¬ 
sence of electromagnetism in Ref. [3] (following Ref. [9]) 
was to apply the twist in a different direction in isospin 
space (ri) from that in which the masses are split (ts). 
This leads to a real quark determinant, and is the method 
used to simulate the s and c quarks using twisted-mass 
fermions (see, e.g.. Ref. uni). This does not, how¬ 
ever, generalize to include electromagnetism in a gauge- 
invariant way. Here, instead, we follow Ref. m, and 
twist in the rs direction. When doing simulations, this 
has the disadvantage of leading to a complex quark deter¬ 
minant]^ but there are no barriers to studying the theory 
with XPT. 

The remainder of this paper is organized as follows. 
We begin in Sec. [TT|with a brief discussion of our power- 


^ These results are from the recent LQCD calculation of Ref. [H], 
and use the convention of that work for the separation of elec¬ 
tromagnetic and Cq effects. 

^ This is avoided in Refs. [Ill [12] by expanding about the theory 
with degenerate quarks and no electromagnetism. 
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counting scheme and a summary of relevant results from 
Ref. [3]. We then explain, in Sec. |III[ how electromag¬ 
netism changes the results of Ref. [J for the case of 
Wilson-like fermions. Section IIVI describes how to si¬ 
multaneously include isospin breaking, electromagnetism 
and twist, while Sec. [V| gives our corresponding results 
for the phase diagram, focusing mainly on the case of 
maximal twist. We conclude in Sec. ED 

Two technical issues are discussed in appendices. The 
first concerns the renormalization factors needed to re¬ 
late lattice masses to the continuum masses that appear 
in XPT. This issue is subtle because singlet and nonsin¬ 
glet masses renormalize differently. This point was not 
discussed in Ref. [3], and we address it in Appendix \A\ 
except that we do not include all the effects introduced 
by electromagnetism. 

The second appendix concerns the need for charge- 
dependent critical masses in the presence of electromag¬ 
netism. These must be determined nonperturbatively, 
and various methods for doing so have been used in the 
literature. One of these methods, proposed in Ref. m, 
can be implemented using partially quenched (PQ) XPT, 
and thus checked. This is done in App. We find that 
the method only provides one constraint on the up and 
down critical masses and must be supplemented by an 
additional condition in order to determine both. 

Appendix requires results from a XPT analysis of 
a theory with twisted nondegenerate charged quarks at 
nonzero lattice spacing and at nonvanishing 6 >qcd- We 
provide such an analysis in a companion paper [13]. 

II. POWER-COUNTING AND SUMMARY OF 
PREVIOUS WORK 

In order to study the low-energy properties of LQCD, 
we must decide on the relative importance of the com¬ 
peting effects. The power counting that we adopt is 


The choice of as the dominant subleading contri¬ 
bution is less obvious, and is discussed in some detail in 
Ref. |3|. The essence of the argument is that, while the 
terms are not necessarily numerically larger than generic 
terms, they give the leading contribution from quark 
mass differences to isospin breaking in the low-energy ef¬ 
fective theory. For example, these contributions give rise 
to the CP-violating phase in the continuum analysis]^ 
In this note we keep only terms up to and including 
those proportional to e^, so that we have the leading order 
term of each type. We refer to this as working at LO+ 
indicating that it goes slightly beyond keeping only LO 
terms. 

We now collect the relevant results from Ref. |3| con¬ 
cerning the phase diagram of Wilson-like fermions in the 
presence of nondegeneracy. We work entirely in SU( 2 ) 
WXPT, in which the chiral field is E G SU(2). The LO+ 
chiral Lagrangian for Wilson-like fermions (whether im¬ 
proved or not) is 

= ( 2 ) 

= -^tr(x^E + Etx) - W'MA^E + Eti)]^ 

+ - Sq)]^ (3) 

where A = 2 IFoal is the spurion field used to introduce 
lattice artifacts. This Lagrangian contains several LECs: 
/ « 92 MeV and Bq from LO continuum XPT, Wq and 
W' introduced by disretization errors, and ij. The latter, 
though of next-to-leading order (NLO) in standard con¬ 
tinuum power-counting, leads to contributions propor¬ 
tional to and thus we keep it in our LO+ calculation. 
£j is not renormalized at one-loop order, and matching 
with SU(3) XPT leads to the estimate [T^ 


m ^ o? q^em > ^ ^ aaEM---, ( 1 ) 

where m represents either rUu or rud. This is the power 
counting adopted in Ref. [3] , except that electromagnetic 
effects are now included. This scheme only makes sense if 
discretization errors linear in a are absent, either because 
the action is improved or because the 0 (a) terms can be 
absorbed into a shift in the quark masses (as is the case 

in WXPT ED- 

The explanation for the choice of leading order (LO) 
terms in this power-counting is as follows. Present sim¬ 
ulations have I/a ~ 3 GeV, and using this together 
with Aqcd ~ 300 MeV we find aAgco ~ O.I. Thus 
second order discretization effects are of relative size 
(aAgco)^ ~ O.OI. This is comparable to q^em: '^u/^qcd 
and rUc^/AgcD (given that rUu ~ 2.5 MeV and rud ~ 
5 MeV [m US]). The results for the neutron-proton 
mass difference described in the Introduction are con¬ 
sistent with this power-counting (using the fact that 
u-rud rud). 


SBoms ’ 

indicating that ^7 is positive. 

The final ingredient in Eq. © is X = 2 ^ 0 M, which 
contains the mass matrix M = diag(m^i, m^^), with 
rUu^d renormalized masses in a mass-independent scheme. 
Since is supposed to represent the long-distance 
physics of a lattice simulation close to the chiral and con¬ 
tinuum limits, to use it we need to know the relationship 
between bare lattice masses and the renormalized masses. 
This relationship is nontrivial when using nondegenerate 
quarks, and is discussed in Appendix This point was 
overlooked in Ref. [3]. 


^ A further justification for this choice, also discussed in Ref. [3], 
is that in SU(3) XPT such terms are of LO, since they are pro¬ 
portional to {mu—md)‘^/ms- 


m. 
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To determine the vacuum of the theory, we must mini¬ 
mize the potential V^. Writing (S) = the potential 

becomes 

= —/^ (mq cos 0 + cie^nl sin^ 0 w' cos^ O) , (5) 

where 


Tflq — 5o (rU-y, , Cq — ‘2tB{)Cq^ 


Q = 


^7 


W = 


MW'W^a^ 


(6) 


Assuming q > 0 [based on the estimate m , the result¬ 
ing phase diagrams are shown in Fig. The unshaded 
phases are continuum-like with |cos^| = 1. The shaded 
(pink) phases violate CP with 


Inal = 1, 


cosO = 



( 7 ) 


The boundaries between continuum-like and CP- 
violating phases he along the lines \rhq\ = 2{c£e^ — w'), 
and are second order transitions. The boundary between 
the two continuum-like phases with opposite cos6> is a 
first order transition. Within the continuum-like phases 
the pion masses are 


mlo = \mq\ - 2{ceel - w'), ml± = \fhq\ + 2w', (8) 

while within the CP-violating phase 

mlo = ‘2{ci'?q - w') sin^ 0 , m^± = 2^?^ . (9) 


The neutral pion mass vanishes along the second or¬ 
der transition lines. Plots of these masses are given in 

Ref. [3]. 


III. CHARGED, NONDEGENERATE WILSON 
QUARKS 

We now add electromagnetism, so that we are consider¬ 
ing Wilson fermions with charged, nondegenerate quarks. 
Precisely how electromagnetism is added at the lattice 
level is not relevant; all we need to know is that electro¬ 
magnetic gauge invariance is maintained by coupling to 
exact vector currents of the lattice theory. We work here 
only at LO in q^em, which in terms of Feynman diagrams 
means keeping only those with a single photon propa¬ 
gator. We also work at infinite volume, thus avoiding 
the complications of power-law volume dependence that 
occur in simulations [iiiTiiiH]. 


fhq 




(b) First-order scenario {w' > 0). 


FIG. 1: Phase diagrams from Ref. [3] including effects 
of both discretization and nondegenerate quarks. CP is 
violated in the (pink) shaded regions. The (blue) lines 
at the boundaries of the shade regions are second-order 
transitions (where the neutral pion mass vanishes), 
while the (yellow) line along the Cq axis joining the two 
shaded regions in ([T^ is a line of first order transitions. 
The analytic expression given for the shaded region in 
(la) holds also for that in (lb). As discussed below in 
Sec. IIIB these phase diagrams apply also in the 
presence of electromagnetism. 


A. Induced shifts in quark masses 


The dominant effect of electromagnetism is a charge 
dependent shift in the critical mass, as noted in Refs. [H 
mi [HI. Here we discuss this shift from the viewpoint 


of the Symanzik low-energy effective Lagrangian [2Ql [21] . 
It arises from QCD self-energy diagrams in which one 
of the gluons is replaced by a photon, and leads to the 
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FIG. 2: Examples of LO contributions from 
electromagnetism to quark self-energies. Diagrams with 
additional gluons and quark loops are not shown. These 
three types of diagram lead, respectively, to the three 
operators listed in Eq. Only the first operator is 
present in the “electroquenched” approximation. 


appearance of the operators 
“ / 

“ /' / 

w^E(4)E7/> (10) 

where f = Cu = 2/3 and = —1/3. Examples of 
the corresponding Eeynman diagrams are shown in Eig.|^ 
These operators are allowed because electromagnetism 
breaks isospin, while Wilson fermions violate chiral sym¬ 
metries. Their contributions are smaller than those of 
the ^fff/a operator that leads to the dominant shift 


in the critical mass. However, because q^em ^ o? ^ m 
in our power-counting, q^em/^^ effects are proportional to 
a ^ and thus dominate over physical quark masses. 

They must therefore be removed by appropriate tuning 
of the bare masses. Since the combined effect of the three 
operators is independent 0 (q^em/<^) shifts in rUu and rrid^ 
removing these shifts requires independent tuning of the 
u and d critical masses. 

Different methods for doing this tuning have been used 
in the literature. The most straightforward, used in 
Ref. [8], is to determine the bare quark masses directly 
by enforcing that an appropriate subset of hadron masses 
agree with their experimental values (keeping all isospin 
breaking effects). This avoids the need to directly deter¬ 
mine the critical masses, but is the most challenging nu¬ 
merically. An alternative approach, proposed in Ref. El, 
makes use of a partially-quenched extension of the theory. 
In Appendix we check this method by showing how it 
can be implemented in XPT. We find that it cannot de¬ 
termine both critical masses, but instead only provides 
a single constraint between them. We then introduce an 
additional tuning criterion which, together with that of 
Ref. m, does allow both critical masses to be deter¬ 
mined. 

Eor the rest of the main text, we assume that the 
charge-dependent critical masses have been determined 
in some manner, such that 0{aEM/ci) self-energy effects 
can be ignored. This leaves electromagnetic corrections 
proportional to q^em, which we must keep in our power 
counting, as well as higher-order effects proportional to 
<^EM X etc., which we can ignore. 

Examples of the latter effects are the bilinears 

<^EM ^ 6/771/// and a^M E ejflPf. (11) 
/ / 

These arise as 0{am) corrections to the operators of 
Eq. ( pR| ), and are also present directly in the continuum 
theory. We stress that, in the Symanzik Lagrangian, one 
has no dimensionful parameters other than m and 1/a, 
so bilinears proportional to q^emAqcd are not allowed. 
Eactors of Aqcd arise when we move from the Symanzik 
Lagrangian to XPT. 

The only effect of electromagnetism that is simply 
proportional to q^em— and thus of LO in our power 
counting—is that arising from one photon exchange be¬ 
tween electromagnetic currents. This is a continuum ef¬ 
fect, long studied in XPT. It leads to he following addi¬ 
tional term in the chiral potential [22j [23] 0 
/2 

Vem = -^cem tr(Er3E’*'T3). (12) 

Here cem is an unknown coefficient proportional to q^em • 
All that is known about cem is that it is positive [24] . 


Contributions from the isoscalar part of the photon coupling lead 
to the same form but with one or both ra’s replaced by identity 
matrices. In either case the contribution reduces to an uninter¬ 
esting constant, and is thus not included in Vem- 
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B. Phase diagram and pion masses 


The competition between electromagnetic effects and 
discretization errors for two degenerate Wilson fermions 
has been previously analyzed in Ref. [25] . Here we add in 
the effects of nondegeneracy. This turns out to be very 
simple. Using the SU(2) identity 


4tr(I]r3Etr3) = [tr(S + St)]" 
together with 


[tr([E-Et]r3T-8, 

(13) 


X = mqH-eqr3, (14) 

we find that Vem can be absorbed into [given in 
Eqs. (§ and <§] by changing the existing constants as 

w' —> w' + cem , and —S- + cem • (15) 

This allows us to determine the phase diagram and pion 
masses directly from the results presented in the previous 
sectionlD 

We first observe that, at the order we work, the phase 
diagram is unehanged by the inelusion of EM —the results 
in Fig. Instill hold. This can be seen from the form of the 
potential in Eq. (§, which, since \ns\ = 1, depends only 
on CiCq — w'. This combination is, however, unaffected 
by the shifts of Eq. ( [T^ and so the phase boundaries 
and values of 0 throughout the phase plane are also un¬ 
changed. 

Similarly, from Eqs. ^ and ^ we see that the neutral 
pion masses are unchanged throughout the phase plane. 
In particular, the second-order phase boundaries are (as 
expected) lines along which the neutral pion is massless. 

The only change caused by electromagnetism is to the 
charged pion masses, which are increased by the same 
amount throughout the phase plane: 

ml± —> ml± + 2cem • (16) 

One implication is that, for = 0, the charged pions are 
no longer massless within the Aoki phase (if present). 
This is because they are no longer Goldstone bosons, 
as the flavor symmetry is explicitly broken by electro¬ 
magnetism. Also, as noted in Ref. [25|, the charged 
pion can be lighter than the neutral one inside the CP- 
violating phases. This is not inconsistent with Witten’s 
identity [24| because the latter did not account for dis¬ 
cretization effects. Plots of the pion masses are shown in 
Fig-lp 

It is perhaps surprising that electromagnetism, which 
contributes at LO in our power-counting, has no effect on 
the phase diagram, whereas the subleading contributions 
proportional to have a significant impact. We can 


^ For 'cq = 0 our results are in complete agreement with those of 
Ref. [25]. 


nK 
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(c) First-order or Aoki scenario with — cem < w' < CiCq 


FIG. 3: Pion masses for nondegenerate untwisted 
Wilson fermions including electromagnetism. The three 
possible behaviors along vertical slices through phase 
diagrams of Fig. are shown. Solid (blue) lines show 
m^o , while dashed (red) lines show . Expressions 
for masses are given in the text. 
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understand this by noting that the CP-violating phase 
is characterized by a neutral pion condensate, which re¬ 
mains uncoupled to the photon until higher order in XPT 
(where form factors enter). 

The implications of these results for practical simula¬ 
tions (such as those of Ref. [8] ) are unchanged from the 
discussion in Ref. [3]. In particular, for the Aoki scenario 
{w' < 0) discretization effects move the CP-violating 
phase closer to the physical point than for degenerate 
quarks, so one must beware of simulating too close to 
this transition. 


IV. NONDEGENERACY, 
ELECTROMAGNETISM AND TWIST 

When using twisted-mass fermions one must decide on 
the relative orientation in isospin space both of the twist 
and the isospin-breaking induced by quark mass differ¬ 
ences and electromagnetism. In the absence of electro¬ 
magnetism, the standard choice is to align these two ef¬ 
fects in orthogonal directions. For example, one usually 
takes T 3 for isospin-breaking, as in the continuum, while 
twisting in the ri directionj^ This is the choice used 
in simulations of the strange-charm sector using twisted- 
mass fermions [26]. It ensures that the fermion determi¬ 
nant is real, and (subject to some conditions) positive |9]. 
This was the choice whose phase structure we determined 
using WXPT in Ref. |3|. 

This approach does not, however, allow for the inclu¬ 
sion of electromagnetism. One problem is apparent al¬ 
ready in the continuum limit, where the twisted-mass 
quark action is (in the “twisted” basis) [27j 

+ rriqC^ -h i^bTirriqS^ -h e^rs)'?/^. (17) 

Here Ip is the gluonic covariant derivative, rriq is the av¬ 
erage quark mass, and uj the twist angle with c^j = cos uj 
and = sin UJ. This action is not invariant under fla¬ 
vor rotations in the rs direction, so there is no conserved 
vector current to which the photon can couple. In other 
words, there is no global flavor transformation available 
to gauge. 

To avoid this problem, we recall that twisting is, in 
the continuum, simply a nonanomalous change of vari¬ 
ables that does not effect physical quantities. Thus we 
should start with the standard action including electro¬ 
magnetism 


'tp{lp - ie4Q + mq + eqTs)^ , (18) 

with the photon field coupling via the charge matrix 

Q = gff + 2^3 5 (19) 


and then perform a chiral twist 

^ ^ ^ ^ . ( 20 ) 

This leads to the quark action of Eq. 0 with the addi¬ 
tion of the photon coupling 


44 


:(Cu;T3 


Sw'r275) 


4- 


( 21 ) 


Thus the photon couples to a linear combination of vector 
currents and to an axial current in the r 2 direction. In 
the continuum, this combination is conserved [given the 
twisted mass matrix of Eq. 0] and the action remains 
gauge invariant. 

We conclude that the correct fermion action to dis¬ 
cretize is the sum of Eqs. 0 and This, however, 

is not possible in a gauge invariant way using Wilson’s 
lattice derivative (except for s^j =0). The Wilson term 
breaks all axial symmetries, so the r 275 part of the pho¬ 
ton coupling is to a lattice current that is not conserved. 

To avoid this problem, and obtain a discretized twisted 
theory that maintains gauge invariance, one needs to 
twist in a direction that leaves the photon coupling to 
a conserved current. The only choice is to twist in the 
Ts direction. Then the twisted form of the continuum 
Lagrangian is 


'lp{lp - ie4Q + rUqCuj + r^CqC^ -h i^^T^rriqS^ -h . 

( 22 ) 

This is discretized by adding the standard Wilson term. 
Since the photon is coupled to vector currents that are 
exact symmetries of both the Wilson term and the full 
mass matrix, gauge invariance is retained. 

This form of the twisted isospin-violating action (with 
UJ = 7r/2) is used in the recent work of Refs. [TTJ [T2] . 
It has one major practical disadvantage—the quark de¬ 
terminant is complex for nonzero twist. This is true for 
nondegenerate masses alone, as explained in Ref. [28] . 
Adding electromagnetism only makes the problem worse, 
since at the least it induces further nondegeneracy in the 
masses. Because the action is complex, direct simula¬ 
tion with present fermion algorithms is challenging. This 
problem is avoided in Refs. [TTl [12] by doing a perturba¬ 
tive expansion in powers of Cq and q^em- The expectation 
values are then evaluated in the theory with no isospin 
breaking, for which the fermion determinant with twist¬ 
ing is real and positive. 

In the following section we study the phase diagram of 
the theory with the discretized form of the Lagrangian 
(22). To our knowledge, this form of the twisted theory 
has not previously been studied in WXPT either with 
nondegeneracy alone or with electromagnetism. 


V. XPT FOR CHARGED, NONDEGENERATE 
QUARKS WITH A T3 TWIST 


® Any linear combination of ri and T 2 is equivalent; ri is the stan¬ 
dard choice. 


The conclusion of the previous section is that the 
twisted-mass theory whose phase diagram is of interest 
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is that with lattice fermion Lagrangian 


the new operators ar^ 


[Dw + mo + Tseo + + ijbVo] • (23) 


is a lattice fermion field and Dw the lattice Dirac 
operator including the Wilson term (and possibly im¬ 
proved). Dw is coupled to both gluons and photons, 
with the latter coupling to the rs vector current. The 
action differs from that considered (implicitly) in Sec. Ill 
only by the addition of the two mass parameters /io and 
^ 0 . 


The four bare mass parameters in (23) are related in 


the continuum limit to the renormalized up and down 
masses, the twist angle (which is a redundant param¬ 
eter) and the QCD theta angle, 6 >qcd- The aim is to 
tune the bare parameters so that the dimension 4 part 
of the quark contribution to the Symanzik effective La¬ 
grangian is given by Eq. (22) with the desired physical 


quark masses, for some choice of uj. As for untwisted Wil¬ 
son fermions the dominant effect of electromagnetism is 
to cause separate 0{aEM/ci) shifts in the (untwisted) up 
and down masses. These shifts depend on twisted masses 
only at quadratic order, so that, to the order we work, 
they are identical to those for Wilson fermions. They can 
be determined by the methods discussed in Sec. |III A] and 
Appendix They are equivalent to independent shifts 
in mo and cq. 

After the additive shift in mo and eo, all four masses 
in (23) must be multiplicatively renormalized in order 


to be related to the continuum masses in Eq. (22). As 


discussed in Appendix this requires different renor¬ 
malization factors for all four masses. We assume here 
that these renormalizations have been carried out, so that 
the dimension four term in the Symanzik effective La¬ 
grangian is given by Eq. (22) and described by the three 
parameters m^, Cq and uj. 

We stress that this tuning and renormalization must 
be carried out with sufficient accuracy. If not, instead of 
Eq. (22), one ends up with a similar form having 


ent twist angles for the rriq and Cq parts. The parity-odd 
parts can then only be removed by a combined flavor 
nonsinglet and flavor singlet twist. Since the latter is 
anomalous, this corresponds to a theory with nondegen¬ 
erate quark masses, electromagnetism, a twist angle and 
a nonvanishing 6 >qcd- In other words, the theory not 
only has the unphysical parity violation due to twisting 
(which can be rotated away in the continuum limit) but 
also the physical parity violation induced by 6 >qcd- In¬ 
deed, to analyze the tuning in XPT one needs to include 
a nonvanishing ^qcd, nn analysis we carry out in a com¬ 
panion paper m- 

Assuming that the dimension-four quark Lagrangian is 
Eq. (22), we next investigate which higher-dimension op¬ 
erators are introduced into the Symanzik Lagrangian by 
twisting. Those operators present for Wilson fermions re¬ 
main, but, as discussed in Sec. |ml are all of higher order 
than we consider. The dominant operators introduced 
by twisting will violate parity, because they are linear in 
the parity-violating mass terms /io and tjq. Examples of 


, aiqo^G^yCr^y^ , and aiiQipr^G^ya^yip . 

(24) 

Since we generically treat am terms as being beyond LO+ 
[see Eq. ([^], we should be able to ignore these opera¬ 
tors. However, because t/o ^ and we are treating Cq as 
somewhat enhanced, one might be concerned about drop¬ 
ping ar/o terms. In fact, the ar/o operators in (24), when 
matched into XPT, pick up an additional factor of m or 
and thus are unambiguously suppressed. The reason 
for the extra factors is that the LO representation of a 
flavor-singlet pseudoscalar in XPT, tr(E — E’*'), vanishes 
identically. Eor the induced ^qcd term, one can also see 
this result by noting that it can be rotated into the isos¬ 
inglet mass term, leading to a contribution proportional 
to mOQCB ^ aem. 

Proceeding in this fashion, we find that all other new 
operators induced by the parity-breaking masses are be¬ 
yond LO+ in our power counting. Thus, once the requi¬ 
site tuning has been done, the LO+ chiral effective theory 
for Ts twisted fermions with isospin breaking is given by 
the same result as for Wilson fermions, i.e. 


f2 

tr + Vx + Vem ( 25 ) 

[see Eqs. ^ and except that the quark mass matrix 
is now twisted 


X = {ffiq + . (26) 

We analyze the phase structure of this chiral theory in 
the next two subsections. 


A. Phase diagram and pion masses at maximal T3 
twist 

We begin working at maximal rs twist, which is the 
choice used in Refs. muig. In this case 


X = ifhqTa + iCq , (27) 

and the chiral potential becomes 

V;^+EM ^ ^ / ^2 I • 2 /) 

-^— = m^ns sm 6 - (c^e^ -\-w )sm 0 

+ cem(cos^ 0 -\-nl sin^ 0 ), (28) 

up to an irrelevant constant. Since cem > 0, the right- 
hand side is maximized always with \ns\ = 1, and we see 
that the cem term becomes a constant. Thus, once again, 
electromagnetism has no impact on the phase diagram. 


^ The first of these corresponds to an induced value of 0 qcd pro¬ 
portional to a?7o • This is one way of seeing that the lattice action 
( |23| > leads to a complex fermion determinant. 
















rriq 



(a) Aoki scenario {w' < 0). 
fhq 



(b) First-order scenario {w' > 0). 

FIG. 4: Phase diagrams including effects of 
discretization and nondegeneracy for maximally 
Ts-twisted quarks. Electromagnetism has no impact on 
the phase diagram. Notation as in Fig. The neutral 
pion is massless along the second-order phase boundary 
between shaded (CP-violating) and unshaded phases. 


We also see that the effect of nondegeneracy can be de¬ 
duced from the results for the degenerate case (studied 
in Refs. [29ff3T] ) simply by shifting w'. 

The resulting phase diagrams are shown in Fig. 
Comparing to the untwisted results of Fig.[^ we see that 
the role of the Aoki and first-order scenarios has inter¬ 
changed. Without loss of generality, we can take = I 
throughout the phase plane. Then, in the continuum- 
like (unshaded) phases we have sinO = sign(mg), cor¬ 
responding to the condensate aligning or antialigning 
with the applied twist. Second order transitions occur 
at \mq\ = 2{w' + C£e^). For smaller values of \fhq\ the 
condensate angle is sin^ = mq/{2[w' + Q?^]), with two 


degenerate minima having opposite signs of cosO. If one 
switches to the “physical basis” in which the twist is put 
on the Wilson term, then one finds that this phase vio¬ 
lates CP, just as in the Wilson case. 

These results differ significantly from the phase struc¬ 
ture for nondegenerate quarks with a maximal ri twist, 
shown in Fig. 3 of Ref. [3]. In particular, an additional 
phase found for w' > 0 with a ri twist is absent here. 
We stress again that only the theory with a rs twist, i.e. 
that discussed here, can incorporate electromagnetism. 

For the pion masses we find the following results. 
Within the continuum-like phases we have 

mlo = \mq\- 2(Qeq + w') , ml± = \rhq\ + 2cem , (29) 

while within the CP-violating phase 

mlo = 2{ceel + w') cos^ 9 , m^± = 2(c<<eq +w' + cem) • 

(30) 

As expected, only the charged pion masses are affected 
by electromagnetism. Plots of these results along vertical 
slices through the phase diagram are shown in Fig. 

It is interesting to compare to the results with un¬ 
twisted fermions, which are given in Eqs. (i and (§ 
together with the shift ( [l^ of by 2 cem induced by 
electromagnetism. We see that the neutral pion mass 
differs only by the change of sign of re' (which also im¬ 
plies the interchange sin^ cos^). This means that 
the results in the two scenarios interchange exactly. For 
the charged pion masses, apart from the interchange of 
scenarios there are also overall shifts proportional to w'. 

The implications of these results for present simula¬ 
tions are as follows. If one could simulate the theory 
directly (somehow dealing with the fact that the action 
is complex) then one would need to avoid working in or 
near the CP-violating phase. This is now more difficult in 
the first-order scenario than the Aoki scenario—opposite 
to the situation with untwisted Wilson fermions. This 
qualitative result is the same as for ri twisting (with¬ 
out electromagnetism), although the area taken up by 
unphysical phases is larger in that case [3]. As noted 
above, actual simulations done to date at maximal twist 
use perturbation theory in Zq and q^em, and so evaluate 
all expectation values in the theory with tq = q^em = 0 - 
Clearly, if re' > 0, these simulations must be careful to 
have mq large enough to avoid the CP-violating phase 


B. Nonmaximal rs twist 

We have also investigated the phase structure for gen¬ 
eral Ts twist, i.e. nonvanishing and nonmaximal. One 


In addition, if these simulations are done close to the onset of 
the CP-violating phase, one would expect the expansion in Pg 
to be poorly convergent. This is probably not a problem for the 
method of Ref. Ea, however, since they take the continuum limit 
of the term linear in 'cq, and in this limit w' = 0 and the lattice 
artifacts discussed here vanish. 
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FIG. 5: Pion masses for nondegenerate maximally 
Ts-twisted fermions including electromagnetism. The 
three possible behaviors along vertical slices through 
phase diagrams of Fig. are shown. Solid (blue) lines 
show m^o, while dashed (red) lines show . 
Expressions for masses are given in the text. 


motivation for doing so is that twisted-mass simulations 
cannot achieve exactly maximal twist; another is to see 
how the phase diagrams of Fig. change into those of 
Fig-i 

Expressions are simplified if we define 0 relative to a 
twist cj, i.e. if we use 


_ ^iujr3/2^ieh-T^iujr3/2 

Then we find (dropping constants)]^ 


(31) 


V 

P 


+ re'(cos ^ cos u; — ns sin^ since)^ 

+ cem(cos^ 6 > + ns sin^ 0) . (32) 


This is not amenable to simple analytic extremization, 
and we have used a mix of analytic and numerical meth¬ 
ods. One can show analytically that the minima always 
occur at |ns| = 1. This implies that, once again, the 
electromagnetism does not play a role in determining the 
phase structure. 

The sign of ns can always be absorbed into 6>, so we can 
again choose ns = 1 without loss of generality. The po¬ 
tential can then be written (up to ^-independent terms) 
as 




= fhq COS 0 + cos^ 0 [—Qe'g -h w' cos(2c(;)] 


713 = 1 


w 


- — sin(2^) sin(2ci;). 


(33) 


A numerical investigation of this potential finds that, for 
nonextremal cj, and for all nonzero re', there is a first- 
order transition as fhq passes through zero, irrespective 
of the value of 7 ^. At this transition 0 jumps from it 12 — 6 
to 7r/2 + (5, with (5^0 depending on the parameters. 
Thus, unlike at the extremal points uj = 0,7r/2, there 
are no second-order transition lines. Correspondingly, 
there are no values of the parameters for which any of 
the pion masses vanish. This is very different from the 
theory with a ri twist, where we found a two-dimensional 
critical sheet [3] in fhq^'eq^uj space. 

The absence of critical lines at nonextremal twist can 
be understood in terms of symmetries. Eor cj = 0 and 
7r/2, the potential has a 6> ^ —6> symmetry, and this 
Z 2 symmetry is broken by the condensate in the CP- 
violating phase, leading to a massless pion at the tran¬ 
sition. Eor nonextremal twist, however, the potential of 
Eq. (32) has no such symmetry. Lacking this symmetry. 


one expects, and finds, only first-order transitions. 


Ata; = 7r/2 this should agree with Eq. ( |28| , and it does once the 
different definitions of 0 and fi are taken into account. 
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VI. CONCLUSIONS 

This work completes our study of how isospin breaking 
impacts the phase structure of Wilson-like and twisted- 
mass fermions. The main results are the phase diagrams 
presented in Figs. [^and[^ together with the correspond¬ 
ing pion masses. These results show how the combina¬ 
tion of discretization errors and nondegeneracy can bring 
unphysical phases closer to (or further away) from the 
physical point. 

The inclusion of electromagnetism into the analysis 
turns out to be very straightforward, aside from the need 
to introduce independent up and down critical masses. 
Electromagnetism has no impact on the phase diagrams 
at leading order, because the condensates in the CP- 
violating phases involve neutral pions. The only impact 
is to uniformly increase the charged pion masses. 

We have investigated within WXPT the conditions 
used in Ref. m to determine the two critical masses in 
the presence of electromagnetism. We find that, unless 
one makes the electroquenched approximation, the two 
conditions are in fact not independent. To determine 
both critical masses one needs an additional condition, 
and we have presented one possibility in Appendixj^ Our 
condition requires simulating at nonzero (though small) 
^QCD, and thus will be difficult to implement in prac¬ 
tice, but provides an existence proof that an alternative 
condition exists. 

Our analysis has been carried out in infinite volume. 
For the finite volumes used in lattice simulations one 
might be concerned about significant finite-volume ef¬ 
fects on the electromagnetic contributions. The impact 
on the results presented here, however, should be mini¬ 
mal. The phase diagram will remain unaffected by elec¬ 
tromagnetism, while the shifts in critical masses are dom¬ 
inated by ultraviolet momenta, themselves insensitive to 
the volume. The only significant effect will be on elec¬ 
tromagnetic mass shifts, with cem picking up an effective 
power-law volume dependence [ilI7l[Il[32]. 
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Appendix A: Relating lattice masses to those in XPT 

In this appendix we describe how bare lattice masses 
used in simulations with Wilson-like fermions are related 
to the masses rUu and rad appearing in XPT (contained 
in the mass matrix M). This discussion draws heavily 
from the results of Ref. [33]. We do not consider the 
impact of electromagnetism here; this is discussed in the 
subsequent appendix. 

We must assume that the number of dynamical quarks 
in the underlying simulations is Nf = 3 (up, down and 


strange) or Nf = 4 (adding charm). Working with up 
and down quarks alone turns out not to be sufficient, but 
in any case this is not the physical theory. We must also 
have that aruf 1 for all flavors /, so that an expansion 
in these quantities makes sense. This condition is met by 
state-of-the-art simulations. Note that this condition is 
much weaker than the requirement that the quarks are 
light in the sense of XPT, which is ruf Aqcd- In the 
main text, we assume the latter condition holds only for 
up and down quarks. 

Let moj be the bare dimensionless lattice mass for 
flavor / (i.e. the mass appearing in the lattice action). 
Because of the additive renormalization induced by ex¬ 
plicit chiral symmetry breaking, unrenormalized quark 
masses are given by 


ruf 


n^o,/ ra^rp 

a 


(Al) 


where mcr is the (dimensionless) critical mass for the 
given number of dynamical flavors. Methods to de¬ 
termine rucr are described below. Then, as shown in 
Ref. [33], renormalized masses are given b}j^ 


ruf 


ffi (^m 1) 


E/ rrif 

Nf 


+ O(am^) 


(A2) 


Here is the renormalization constant for flavor non¬ 
singlet mass combinations such as Cq = {rUu — rud) 
while ZjnTm is the corresponding constant for the aver¬ 
age quark mass, — 1 is a finite constant, arising first at 
0(^^) in perturbation theory. By implementing contin¬ 
uum Ward-Takahashi identities, one can determine Vm 
nonperturbatively for Nf = 3 and 4, although not for 
Nf = 2 [33]. This is the reason for the restriction on Nf 
noted above. We assume here that Vm has been calcu¬ 
lated in this way . 

Equation (A2) shows that the renormalized mass ruf 
does not vanish when rfif = 0 if other flavors are massive. 
Specifically, for the up and down quarks we have 


mu^md = Zm —^— [mu^md) 

+ Zm {fhs + fhch ), (A3) 

mu-md = Zm{mu - md) • (A4) 


(Here we we have chosen Nf = 4 for definiteness; the re¬ 
sult for Nf = 3 is similar.) Thus the two-flavor massless 
point receives an overall additive shift due to the strange 
and charm quarks, and we also see explicitly the differ¬ 
ence between singlet and nonsinglet renormalizations. 

These results imply that, in terms of unrenormalized 
masses, the phase diagrams of Fig. would be trans¬ 
lated in the vertical direction (due to the additive mass 


The correction terms of 0{arh^) in ( |A2| ) are subleading in our 
power-counting and will be dropped henceforth. 
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shift) and stretched by different factors in the vertical and 
horizontal directions. The respective stretch factors are 
+ ^m)/2 and If, however, is known, 

then the two stretch factors can be made equal by ap¬ 
plying a finite renormalization to remove the (1 + r^)/2 
factor. Knowledge of is, however, not useful, since it 
always appears multiplied by the unknown LEG Bq. 


We would like to be able to remove the additive mass 
shift in Eq. (A3). To do so we consider how the critical 
mass TUcr is determined. The expressions above assume 
that it has been obtained by doing simulations with Nf 
degenerate quarks of mass m, and equating nicr to the 
value of m at which the “PCAC mass” vanishes. This is 
equivalent to imposing 


{Tr+\df,{u-t^j5d)\0) 


= 0 . 

m=mcr 


(A5) 


If, instead, one imposes this condition by varying m = 
rUu = with rrig and rrich held fixed at their physical 
values, then the rucr so obtained automatically includes 
the shift due to loops of strange and charm quarks. This 
is because one is enforcing a consequence of chiral sym¬ 
metry in the two-flavor subsector. With this new choice 
of rUcr, and with the adjustment of stretch factors de¬ 
scribed above, the phase diagrams of Eig.[2apply directly 
for lattice masses frif. 

This new choice of nicr has a second advantage: it re¬ 
moves an additional shift of 0(a) in the relation between 
bare quark masses and the masses appearing in XPT. 
As explained in Ref. [2], this shift is caused by the 0{a) 
clover term in the Symanzik effective action (and is thus 
absent for nonperturbatively improved Wilson fermions). 
In the main text it is assumed that this shift has been 
removed. 

Since we include O(a^) terms in the main text, we must 
determine how they impact the considerations above. 
There is no further shift in the quark masses at this 
order—this next occurs at 0{a^) [34]. However, as il¬ 
lustrated by Eig.[^ the O(a^) terms do impact the phase 
diagram. This means that, in general, one cannot use 
the vanishing of the PC AC mass to determine rUcr with 
untwisted Wilson fermions. Eor example, if one is in the 
first-order scenario [Eig. I^b) along the ifiq axis], then 
the PC AC mass simply does not vanish for any In¬ 
stead, one must introduce a twisted component to the 
mass, /i ^ 0{a), and then enforce the vanishing of the 
PC AC mass (in the so-called “twisted basis”). Extrapo¬ 
lating the result linearly to /i = 0 yields a result for rUcr 
that has errors of O(a^), which is sufficiently accurate for 
our analysis. Eor a detailed discussion of this point see 
Ref. [34|. 

In summary, by determining from Ward identities, 
and the critical mass from the PC AC mass condition with 
twisted-mass quarks, one can obtain lattice quark masses 
which are proportional to those appearing in XPT at the 


order we work. Specifically, we find 


rua 




BoZrr 


-{rhu+'md) and 


= imu-fhd), 

noZm 

(A6) 

where rhq and are the quantities appearing in the chiral 
potential of Eq. (§. 

This analysis can be straightforwardly extended to ar¬ 
bitrary twist. We begin with maximal twist, for which 
the mass matrix in XPT is given by Eq. (27), and the 
relevant bare masses are /io and r]o of Eq. (23). In this 


case there is no additive renormalization, but the pres¬ 
ence of different renormalization factors for singlet and 
nonsinglet masses remains. Using the results of Ref. [33] . 
we fin(£3 


rUn 


BnZ. 


O^m 


Zs l+rp Cq Zs ... 

^ ;—Aio and —— = —% . (A7) 

Zp rp 


Here ZsjZp and rp are finite constants, both of which 
can be determined from Ward identities for Nf = 3 and 
4, but not for Nf = 2 [33|. Like r^^ 'f'p begins at 0{g^) 
in perturbation theory. 

At arbitrary twist one has four bare masses, and 
they are related to the corresponding four renormalized 
masses using the same renormalization factors as given 
in Eqs. (A6) and (A7). 


Einally, we stress that the analysis presented here does 
not include electromagnetic effects. The dominant such 
effect is that the critical mass rUcr has to be chosen differ¬ 
ently for the up and down quarks, and is discussed in the 
following appendix. A subdominant, but still important, 
effect is that the renormalization factors now depend not 
only on as but also on q^em- The latter dependence can 
presumably be adequately captured using perturbation 
theory. The formulae given above still hold if one uses 
the new critical masses and renormalization factors. 


Appendix B: Determining the critical masses in the 
presence of electromagnetism 

The analysis of the previous appendix must be ex¬ 
tended when electromagnetism is included, due to the 
presence of charge-dependent self energy corrections pro¬ 
portional to aEu/d- This implies that the critical masses 
for up and down quarks differ, and we label them mcr,u 
and TUcr4^ respectively. In Ref. m two methods for a 
nonperturbative determination of these critical masses 
are proposed. One of these (the method used in prac¬ 
tice in Ref. m) involves only up and down quarks, and 
thus can be implemented, and therefore checked, within 
SU(2) WXPT. We do so in this appendix, finding that the 
method does not fix both eritieal masses, but rather con¬ 
strains them to lie in a one-dimensional subspace of the 


Specifically, we have used Zm = ^lZs and Vm = 1/^S'- 
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^cr,u— '^cr.d plane. We then provide an additional con¬ 
dition that does completely determine mcr,u and rricr.d- 
The tuning conditions require the use of twisted-mass 
quarks, although the resulting values of rricr^u and mcr,d 
apply for both Wilson and twisted-mass quarks. Thus 
the lattice quark Lagrangian is given by Eq. ( [^ . We 
can write the mass matrix in two useful forms 


mo + Tseo + ijsrsido + 

m^OjU T 0 

0 mo,d - 

The tuning proceeds by first choosing bare twisted 
masses /io,u and /io,d such that, when multiplicatively 
renormalized as described in the previous appendix, they 
give rise, respectively, to the desired physical up and 
down quark massesThe negative sign multiplying /io,d 
is chosen to correspond to a rs twist. The second step 
is to tune the untwisted masses mo,u and mo,d to their 
critical values such that the (additively) renormalized un¬ 
twisted masses vanish. 

The method of determining nicr used in the previous 
section is no longer useful—the vanishing of the PC AC 
mass is a condition based on the recovery of the chiral 
SU(2) group, but this group is explicitly broken by elec¬ 
tromagnetism. The workaround proposed in Ref. m 
is to add to the sea quarks (labeled us and ds) a pair 
of valence quarks, uy and dy^ each of which has the 
same charge and untwisted mass as the corresponding sea 
quark, but has opposite twisted massp] Thus {us^uy) 
and {dy^ds) each form a twisted pair. The key point 
is that, within each pair, the 0(<aEM/u) shift in the un¬ 
twisted mass is common. Therefore it is plausible that 
one can determine the critical mass for each pair by en¬ 
forcing the recovery of the corresponding valence-sea chi¬ 
ral SU(2). Specifically, mcr,u is determined by 




{T^sv\9,,{usl,,'y5Uv)\0) 


= 0 , 

'^0, u — 'TTlcr, u 


(B2) 


while mcr,d is determined by the analogous condition 
with u ^ d: 


{T^sv\dnids'lf^'l5dv)\0) 


= 0 . 

'^0 ,d—'^cr,d 


(B3) 


Here ir'^y and are sea-valence pions composed, re¬ 
spectively, of up and down quarks. 


In fact, the tuning can be done using any values of the twisted 
masses which respect our power counting. The critical masses 
do not depend on the twisted masses at the order we work. 

This description is equivalent to that of Ref. m, but differs tech¬ 
nically in two ways. First, we find that one need only introduce 
two valence quarks to describe the method, rather than the four 
used in Ref. m- This does not impact the method itself, only 
its description. Second, we work in the twisted basis, rather than 
the physical basis used in Ref. HU. 


When using a partially quenched theory, one also needs 
to add ghost fields, uy and dy, to cancel the valence 
quark determinant sp^ Thus the full softly-broken chi¬ 
ral symmetry is the graded group SU{4\2)l x SU{4\2)r. 
This raises the question of whether complications arising 
from partial quenching, or from discretization eff ects , can 
lead to corrections to the tuning criteria of Eqs. ( |B2[ ) and 
(B3). This is one of the issues we address here by map¬ 


ping these conditions into XPT. 

We begin by mapping the mass matrix in the un- 
quen ched sector into XPT. The four parameters of 
Eq. (Bl) map into 


X = 


rriue 

0 


0 

-iujd 


rride 

0 


( to , - 


(B4) 

(B5) 


The choice of sign for ujd is such that it is positive with a 
Ts twist. X contains the additional parameter ip compared 
to the mass matrix analyzed in the main text, Eq. (26). 


(/? is a measure of the difference between up and down 
twist angles. 


Ldu = ^ ^, ujd = ^ ^ • (B6) 


As discussed in Sec. |V| such a difference corresponds to 
the introduction of a nonzero Oqcb —the explicit relation 
is if = 6>qcd/2. 


We no te that the relations between th e ba re masses of 
Eq. (Bl) and the parameters of y in Eq. (B5)—which can 
be worked out along the lines of the previous appendix— 
are not needed here. All we need to know is that, if 
mo^u = mcr,u and mo,d = mcr,d^ then both up and down 
masses are fully twisted. Thus the twist angles in y are 
= ^d = 7r/2, implying maximal twist with no ^qcd 
term: uj = Trj2 and (p = 0. Reaching this point in pa¬ 
rameter space is the aim of tuning. 

When considering the PQ extension of this theory, we 
will focus mainly on the quark sector, since the ghosts 
do not play a significant role. Collecting the four quark 
fields in the order 


'^PQ = uv, dy, ds ), (B7) 

the extended quark mass matrix is 


Xpq = 


{fhq -h 

0 {fhq 


eg)e^‘ 


(B8) 


The factors of rs arise because, by construction, valence 
quarks have opposite twisted masses to the correspond¬ 
ing sea quarks. We stress that the 0{ay,M/o) shifts are 
incorporated into the parameters fhq and Tg, along with 


16 Pqj. i-eviews of partially quenched theories and the corresponding 
XPT, see Refs. [35] [36]. 
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the usual 0{l/a) shifts. We can also include the 0{a) 
shifts in the same fashion. 


To implement the conditions (B2) and (B3) in the PQ 
theory, we need the extension of S to this theory. This 
is a 6 X 6 matrix transforming in the usual way under 
SU (4|2) L X SU (4|2)i^. In fact, as we only need matrix el¬ 
ements for states composed of quarks, and since we know 
from Ref. [37] that there are no quark-ghost condensates, 
we can focus on the 4x4 quark sub-block, which we call 
SpQ- We now argue that the expectation value of T^pq 
has the form 


(Spq) = diag(e®'',e . 


(B9) 


This is based on the following results. First, the un¬ 
quenched 2x2 block of SpQ (i.e. that involving the 
first and last rows and columns) is just the unquenched 
S field. This is unaffected by partial quenching [38l [39] . 
and its expectation value is given by an unquenched XPT 
calculation. This calculation must include not only non¬ 
degeneracy, electromagnetism and twist, but also nonva¬ 
nishing 6>qcd- To our knowledge such an analysis has not 
previously been done, so we carry it out in a companion 
paper m- The result is that the unquenched conden¬ 
sate (H) only rotates in the ts direction—there are no 
off-diagonal condensates suc h as {usds)- This fixes the 
first and last entries in Eq. ( |B9| ) to have opposite phase 
angles. 

This is an important result for the following, so we 
emphasize its key features. Although ^qcd 7^ 0 leads to 
an overall phase in the mass matrix in Eq. (B5)], its 
effect on the condensate (E) is qualitatively similar to 
that of a twist cj, despite the fact that the latter leads to 
opposite phases on u and d quarks. This happens because 
E is constrained to lie in S'1/(2), and so has no way to 
break parity other than rotating in the ts direction. An 
overall phase rotation would take it out of SU (2) into the 
1/(2) manifold. 

The second result needed to obtain Eq. ( |B9| ) is the ex¬ 
istence of relations between valence and sea-quark con¬ 
densates. In particular, one can show that 


more valence quarks, e.g. {uydy) and {uyds)- These dif¬ 
fer from the diagonal condensates in that there is no mass 
term in the quark-level Lagrangian that can serve as a 
source for such condensates. Thus to determine whether 
they are nonzero one must add a source, e.g. AdyUv^ 
calculate the resulting condensate, send the volume to 
infinity, and finally send the parameter A ^ 0. This 
analysis has been carried out in Appendix A of Ref. HOI 
in a theory with twisted-mass quarks, although, unlike 
our situation, the quarks were degenerate and Oqcb = 0. 
The general lessons from Ref. [40] are (i) that to ob¬ 
tain a nonvanishing condensate one needs a source of 
infrared divergence to cancel the overall factor of A, and 
(ii) that nonvanishing twisted masses cut off such diver¬ 
gences. These lessons apply also for all the off-diagonal 
condensates that we consider here. However, the argu¬ 
ment as given in Ref. [40] assumes that the measure is 
real and positive, which does not hold here. Neverthe¬ 
less, since we are tuning to ^qcd = 0 , we expect the 
impact of having a complex measure to be small. Eur- 
t her more, we know from Ref. m that the correspond¬ 
ing sea quark condensates, e.g. {usds) and (us^^ds), 
vanish even when ^qcd 7 ^ 0. These condensates differ 
from those containing valence quarks only by changing 
the signs of some of the twisted masses. Since it is the 
presence of these masses, and not their detailed proper¬ 
ties, that leads to the vanishing of the condensate, we 
expect the result hold s for all off-diagonal condensates. 

With the form (B9) in hand, we can now apply the tun¬ 
ing conditions (B2) and (B3) in XPT. We do so by gen¬ 


eralizing the analysis of Ref. m, where the twist angle 
for unquenched twisted-mass fermions was determined in 
XPT by applying a PC AC-like condition. The required 
extension is from the SU (2) sea-quark sector alone to the 
full valence-sea 5'//(4) symmetry. Much of the analysis 
carries over with minimal changes from Ref. m , so we 
only sketch the calculation. 

The first step is to obtain the pion fields that couple 
to external particles in the tuning conditions. Eollowing 
Ref. [41], we obtain these by expanding the chiral field 
about its vacuum value as 


(uyuv) = (usus) and {uv75Uv) = -{uslbUs), (BIO) 


to all orders in the hopping parameter expansion. The 
minus sign in the second relation follows from the op¬ 
posite twist ed mass of sea and valence quarks. The 
result (BIO) holds on each configuration and thus also 


for the ensemble average, even though the measure is 
complex for ^qcd 7 ^ 0. Since the additive and multi¬ 
plicative renormalizations of these conden sates are the 
same for valence and sea quarks, the result (BIO) implies 


that valence and sea up-quark condensates have opposite 
“twists”, The same argument applies to the down- 

quark condensates, and taken together these arguments 
determine the form of the second and third diagonal el¬ 
ements in Eq. (B9). 


The final result needed to obtain the form (B9) is the 


vanishing of off-diagonal condensates involving one or 


'PQ = ipQ , 

(Bll) 

15 


n = y] TTaAa , 

a=l 

(B12) 


^PQ = = diag(e'^/^, e e *^/^). 

(BI3) 


Here are the generators of SU(4), with tt^ the cor¬ 
responding pion fields. These are the pions in the PQ 
theory that are composed of quarks alone, with no ghost 
component!^ The pions needed for tuning, p'gy and 


A similar form to Eq. ( |B13| > holds for the full 6x6 PQ chiral field, 
but we can focus on the SU{A) block, since the pions we leave 
out in this way are those containing one or more ghost fields. 
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TT^y, are contained in the upper and lower diagonal 2x2 
blocks of n, respectively. 

The next step is to determine the form, in XPT, of the 
axial currents appearing in the tuning conditions. These 
can be obtained by introducing sources into derivatives 
using standard methodology. Since, by definition, our 
chiral potential does not contain derivatives, at LO+ only 
the LO kinetic term [shown in Eq. (§] enters into the 
determination of the currents. We do not display the 
form of the currents, however, as the calculation needed 
for each of the tuning conditions is exactly the same as 
that carried out in Ref. m- This is because each tun¬ 
ing condition involves a separate, nonoverlapping SU{2) 
subgroup of SU (4) (upper-le ft or lower-right 2x2 block), 
and because the condensate (B9) does not connect these 
subgroups. We simply quote the results of the calcula¬ 
tion: 


{t^ sv\^^l{usiiil^uv)\Q) cos 0, (B14) 

sv\dn{dsini^dv)\0) (X cose. (B15) 

Thus enforcing either ( |B2[ ) or ( |B3| ) has the effect of set- 
ting 0 = ±7r/2 and the condensate to 

(Spg) = ±diag(i, -i, i, -i ), (B16) 


For our choices of signs of the twisted masses /io,w and 
/io,d in Eq. (Bl), the ± signs are in fact plusses, i.e. 0 = 
7r/2. 

A surprising aspect of this result is that the two tun¬ 
ing c ond itions are n ot in dependent: if one enforces, say, 
Eq. (B2) then Eq. (B3) will be automatically satisfied. 
This dependence arises because changing rUu in turn 
changes ip and uj and this impacts the d condensate 
through the quark determinant. One might, therefore, 
wonder how the two tuning conditions have been suc¬ 
cessfully applied in Ref. m- To understand this, we note 
that this work makes two approximations. First, isospin¬ 
breaking effects are evaluated only through linear order 
in an expansion in rUu—md and q^em- Second, insertions 
of rriu — md or photons on sea-quark loops are dropped 
(the “electroquenched approximation”). The latter ap¬ 
proximation has the effect of disconnecting the two tun¬ 
ing conditions—all quark loops in both conditions are 
evaluated with uncharged, degenerate sea-quarks, so the 
rt-quark condensate cannot be impacted by changes in rad 
and vice versa. Since XPT predicts that there is a tight 
correlation between the condensates, it appears to us that 
the electroquenched approximation is theoretically prob¬ 
lematic. However, from a purely numerical viewpoint, 
the dropped contributions may well lead only to small 
corrections. 


The lack of independence implies that the tuning con¬ 
ditions cannot determine both rUu^c and md,c—only one 
constraint on these two critical masses is o btained. In 
terms of the parameters of mass matrix (B5), the condi¬ 
tions determine only a rela tion bet ween uj and p. Thus, 
after enforcing either (B2) or (B3) the theory is known 


to lie along a line in the uj— p plane. In terms of the 
bare masses, the theory lies along a line in the mo,u— 
niQ^d plane (with, recall, /io,w and /io,d fixed at the values 
leading to physical quark masses when mo,u = = 0)- 

We do know that this one-dimensional subspace includes 
the point we are trying to tune to, namely that with 
(uj^p) = (7r/2,0). This follows from the analysis of 
Sec. |V A[ At maximal twist with p = 0^ the twist in 
the condensate is also maximal, i.e. 0 = 7r/2. The only 
caveat is that the values of the twisted masses must be 
such that one lies in the continuum-like phase, rather 
than the CP-violating phase (see Fig. [^. 

To complete the tuning we need an additional condi¬ 
tion that forces us to the desired point along the allowed 
line. At first blush one might expect that it would be 
simple to find an additional tuning condition, since theo¬ 
ries with Oqcb 7^ 0 have explicit parity violation. This is 
in contrast to the parity violation induced by a nonzero 
twist UJ which, in the continuum limit, can be removed 
by a chiral rotation. This suggests that one should look 
for quantities that vanish when parity is a good symme¬ 
try. The flaw in this approach is that parity is broken by 
cj 7 ^ 0 away from the continuum limit—the chiral rota¬ 
tion required to obtain the parity-symmetric form is not 
a symmetry on the lattice. Thus the distinction between 
p ^ 0 and cj 7 ^ 0 no longer holds. 

The only choice that we have found for a second con¬ 
dition involves using the pion masses. Specifically, we 
find that, along the line picked out by setting 0 = 7r/2, 
the masses of both charged and neutral pious are mini¬ 
mized when p = 0. This assumes only that we are in the 
continuum-like phase for the physical values of /io,u and 
d'0,d- 

The details of the calculation are presented in Ref. [13] . 
Working at LO+, we find that the constraint 0 = nj2 
forces the quark masses to lie on the line 

rhu + ceifiu - M J ' 


As noted above, this line passes through the desired point 
rUu = rad = 0- The slope is —1 when q = 0, and in¬ 
creases in magnitude as q increases (assuming the phys¬ 
ical situation pu < Pd)- There is no singularity when the 
slope reaches infinity—this simply means that the con¬ 
straint line is the rau = 0 axis. For larger q the slope is 
positive. It decreases with increasing q, though it always 
remains greater than unity. The pion masses along the 
constraint line are 


mlo 


ml± 


f^u T pjd 


- 2q 


t^d 


+ 2q 


mu-rad 


+ 2q 


mu-md 


-2w', 
2 


+ 2cem • 


(B18) 

(B19) 


Thus we see that both masses are minimized along the 
constraint line when mu = md = 0. If one were to imple- 
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ment this tuning condition in practice, then one would 
apply it for the charged pion masses, since these have no 
quark-disconnected contractions. 

This analysis breaks down when q gets too large, be¬ 
cause the theory with rUu = rrid = 0 then lies in the 
the C P-violating phase. This can be seen from the result 
(B18)—for large enough q the squared neutral pion mass 


if one perfectly implements our two tuning conditions— 


becomes negative. This happens sooner for the first-order 
scenario, w' > 0. 


namely either Eq. (B2) or (B3) and minimizing the pion 
masses—one will not have precisely tuned to rriu = = 

0. This can be seen, for example, from the analysis of 
Ref. [41], where terms of 0{ap^,am) lead to maximal 
twist occurring at untwisted masses of 0{aii)^ with /r 
the twisted mass, rather than zero. Shifts of this size 
occur also in the presence of isospin breaking, although 
the detailed form of the corrections will differ. Within 
our power-counting, however, /i ^ so that the shifts 


We close this section by commenting on the impact of 
higher-order terms in APT. Because of such terms, even 

in the untwisted masses are ^ beyond the order 

that we consider. 
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